On a common circle: natural scenes and Gestalt rules 
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Abstract 

To understand how the human visual system analyzes images, 
it is essential to know the structure of the visual environ- 
ment. In particular, natural images display consistent statis- 
tical properties that distinguish them from random luminance 
distributions. We have studied the geometric regularities of 
oriented elements (edges or line segments) present in an ensem- 
ble of visual scenes, asking how much information the presence 
of a segment in a particular location of the visual scene car- 
ries about the presence of a second segment at different rela- 
tive positions and orientations. We observed strong long-range 
correlations in the distribution of oriented segments that ex- 
tend over the whole visual field. We further show that a very 
simple geometric rule, cocircularity , predicts the arrangement 
of segments in natural scenes, and that different geometrical 
arrangements show relevant differences in their scaling prop- 
erties. Our results show similarities to geometric features of 
previous physiological and psychophysical studies. We discuss 
the implications of these findings for theories of early vision. 

One of the most difficult problems that the visual 
system has to solve is to group different elements of 
a scene into individual objects. Despite its computa- 
tional complexity, this process is normally effortless, 
spontaneous and unambiguous 0. The phenomenol- 
ogy of grouping was described by the Gestalt psy- 
chologists in a series of rules summarized in the idea 
of good continuation ||, More quantitative psy- 
chophysical measurements have shown the existence 
of association fields [01 or rules that determine the 



interaction between neighboring oriented elements in 
the visual scene JH[ |(|. Based on these rules and 
on the Gestalt ideas, pairs of oriented elements that 
are placed in space in such a way that they extend 
on a smooth contour joining them will normally be 
grouped together. 

These psychophysical ideas have been steadily 
gaining solid neurophysiological support. Neurons 
in primary visual cortex (VI) respond when a bar 
is presented at a particular location and at a spe- 
cific orientation j?J. In addition, the responses of 
VI neurons are modulated by contextual interactions 
§, H |, 0, 0, |l| 0, |l|, such as the joint 
presence of contour elements within the receptive 
field and in its surround. This modulation depends 
upon the precise geometrical arrangement of linear 
elements [l6), in a manner corresponding to the 
specificity of linkage of cortical columns by long-range 
horizontal connections \ V7, |l^|. Thus neurons in VI 



interact with one another in geometrically meaning- 
ful ways, and through these interactions, neuronal 
responses become selective for combinations of stim- 
ulus features that can extend far from the receptive 
field core. 

The rules of good continuation, the association 
field, and the connections in primary visual cortex 
provide evidence of interaction of pairs of oriented 
elements at the psychophysical, physiological and 
anatomical level. The nature of the interaction is de- 
termined by the geometry of the arrangement, includ- 
ing spatial arrangement and the orientation of seg- 
ments within the visual scene. An important question 
is whether this geometry is related to natural geomet- 
ric regularities present in the environment. It is well 
known that natural images differ from random lumi- 
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nance distributions p0[ but the structural studies 
of natural scenes have not yet addressed the existence 
of geometrical regularities. Here, we address this is- 
sue, by studying whether particular pairs of oriented 
elements are likely to co-occur in natural scenes as 
a function of their orientation and relative location 
in space. Our results are focused on two different 
aspects of the organization of oriented elements in 
natural scenes: scaling and geometric relationships. 
We will show that these two are interdependent. 

Scaling measurements involve studying how the 
probability of finding a co-occurring pair changes as 
a function of the relative distance. A classic result in 
the analysis of natural scenes is that the luminance of 
pair of pixels is correlated and that this correlation is 
scale invariant (l9[ |2pj . This indicates that statisti- 
cal dependencies between pair of pixels do not depend 
on whether the observer zooms in on a small window 
or zooms out to a broad vista. The scale invariance 
results from stable physical properties such as a com- 
mon source of illumination and the existence of ob- 
jects of different sizes and similar reflectance proper- 
ties Here we show that for particular geometries, 
the probability of finding a pair of segments follows 
a power law relation and thus is scale invariant. We 
show further that a very simple geometric rule, con- 
sistent with the idea of good continuation, predicts 
the arrangement of segments in natural scenes. 

Materials and Methods 

Images were obtained from a publicly available 
database (http ://hlab .phys. rug .nl/imlib /index. html) 
|p2| of about 4000 uncompressed black and white pic- 
tures, 1536x1024 pixels in size and 12 bits in depth, 
with an angular resolution of "1 minute of arc per 
pixel. This particular database was chosen because of 
the high quality of its pictures, especially the lack of 
motion and compression artifacts, which would oth- 
erwise overwhelm our statistics. To obtain a mea- 
sure of local orientation we used the steerable fil- 
ters of the H 2 and G*2 basis [||). Using steerable 
filters, the energy value at any orientation can be 
calculated by extrapolating the responses of a set 
of basis filters. A G2 filter is a second derivative 
of a Gaussian and the H2 filter is its Hilbert trans- 



form. H2 and G2 filters have the same amplitude 
spectra, but they are 90° out of phase; that makes 
them a quadrature pair basis filters. The size of the 
filters used was 7x7 pixels. A measure of oriented 
energy was obtained by combining both sets of filters 
E(tp) = Gl(ip)+H2((p) Q. This measure is repeated 
at every pixel of the image to obtain the energy func- 
tion for each image (n) of the ensemble {E n (x, y, (p)}. 
To study the joint statistics of E(x,y,tp), we dis- 
cretized the different orientations at 16 different val- 

np« D = ( =Z- JL\ 1 = (JL 3JL\ 1^ = (2^JL 31lL) 

uca, kj V 32 ' 32'' ^32' 32'' '"' \ 32 ' 32 / 

as shown in the color representation of orientations 
of Fig. 1. With this information one can obtain a 
measure of the statistics of pairs of segments by cal- 
culating the correlation (weighting the co-occurrences 
of segments by their energy) 

C{Ax,Ay,y,ip) = 



Jj^2 E n (x,y,ip)E n (x + Ax,y + Ay,if>)dxdy 

n=l 

where TV is the total number of images and the inte- 
gral is over each of the images of the ensemble. We 
were interested in measuring long-range correlations 
so we studied values of Ax, Ay = {—256,256}. The 
correlation matrix has dimensions 512x512x16x16 
and each point results from averaging 4000 integrals 
over a 1536x1024 domain. To simplify the compu- 
tations, for the general case, we decided to store 
at each pixel, for every image, the maximum en- 
ergy value E(ip max ) and its corresponding orienta- 
tion ifmax- An energy threshold Et was arbitrar- 
ily set to match the visual perception of edges in 
a few images. Pixels in an image were considered 
"oriented" if E({p max ) > Et, "non-oriented" other- 
wise. This unique threshold value was applied to 
all images in the ensemble. Thus, for each im- 
age, we extracted a binary field E^ n (x,y) = {0,1} 
and an orientation field Ang n (x,y) = {1, 16}. 
From this binary field we can construct a histogram 
of co-occurrences: how many times an element at 
postion (x,y) was considered oriented with orien- 
tation if and at position (x + Ax,y + Ay) a seg- 
ment was considered oriented with orientation ip. 
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Thus, formally, the histogram is obtained as C, tak- 
ing as the Energy function E n (x,y,ip) = 1 if ip = 
Ang n (x,y) and E% n (x,y) = 1; E n (x,y) = in any 
other case. The computation is reduced to counting 
the co-occurrences in the histogram ff (Ax, Ay, ip, ip) 
with Ax = {= 256,256},Ay = {-256, 256},^, V = 
(0, tt/16, 27r/16, 7r). From the histogram we ob- 
tained a measure of statistical dependence. While 
choosing the threshold followed computational rea- 
sons, cortical neurons perform a thresholding op- 
eration and thus the measure of linear correlation 
(weighting co-occurrences by their energy) is not nec- 
essarily a more accurate measure of statistical depen- 
dence. The histogram was used for all the data shown 
in Fig. 2a, 2b and 2c, Fig. 3, Fig. 4, and Fig. 5. For 
Fig. 2d, for the particular case of collinear interac- 
tions we computed the full linear cross-correlation. 
This computation is considerably easier since it is 
done for fixed values of orientation and direction in 
space. The two measures shown (Laplacian correla- 
tion and collinear correlation) were obtained accord- 
ing to the formulas: 

C(r) = 

J2 J2 E Lap {x iy )E Lap (x + Ax,y + Ay) 

x,V Ax 2 +Ay 2 =r 2 

-(j2E Lap (x, y )r 

for Laplacian filtering, and 

C{r) = ]T E(x, y, 0)E(x + r, y, 0) - (£, E ( x > ^ °)) 2 

x,y x,y 

for collinear oriented filtering. 

A quantitative signature of scale invariance is given 
by a function of the form C = r~ a (power law) where 
C is the correlation, r the distance and a constant. If 
the scale is changed r — > Ar = r' the function changes 
as C(r) = \~ a r~ a = fcC(r') where A; is a constant. 
A power law is easily identified as a linear plot in 
the log-log graph, which is clear from the relation 
log(C) = -alog(r). 



The axis of maximal correlation (Fig. 5b) was 
calculated as follows. For each pair of orientations 
((/?, ip), a measure of co-occurrence was calculated in- 
tegrating across 16 different lines of angles of val- 
ues (0, jfg, — i i") over distances of [-40,40] of the 
center of the histogram. Thus, for an angle 9 and 
orientations (ip, ip) the measure of co-occurrence is: 

P<pA b ) = EiL 4 o#( cos (0) * * i,tp,1>)- We 

then calculated the direction of maximal correla- 
tion max (p,ip) and we grouped all angles with com- 
mon relative orientation, tp — ip = £. We had 16 
different values for each £ and from these 16 dif- 
ferent values we calculated the mean P(0, e) =< 
9max{4'j' l P + e) >y> and the standard error. To calcu- 
late the mean energy as a function of relative orien- 
tation (Fig. 3) we integrated the histogram in spa- 
tial coordinates for each pair of orientations in space, 
and, as before, the different pairs where grouped 
according to their relative difference in orientation 
to calculate a mean value and a standard devia- 
tion, E v ^ = J x= _ 1QQ f y= _ WQ H(x,y,(p,ip)dxdy and 
E(ip) = (E e ^ +lp ) e . The code was parallelized using 
MPI libraries and run over a small Beowulf cluster of 
Linux workstations. 

In general, horizontal and vertical directions had 
better statistics since there are more horizontal or 
vertical segments than oblique in the images; these 
special orientations are also the ones most prone to 
artifacts from aliasing, staircasing, and the ensem- 
ble choice. Since we are interested in this study in 
the correlations as a function of relative distance and 
orientations, all the quantitative measurments were 
performed averaging over all orientations. However, 
the results shown still held true for each individual 
orientation. 

Results 

All 4000 images used in this study were black and 
white, 1536x1024 pixels in size and 12 bits in depth. 
We used a set of filters to obtain a measure of orienta- 
tion at each pixel of every image of the database |p3| . 
The filters were 7x7 pixels in size and thus provided a 
local measure of orientation. The output of the filter 
was high at pixels where contrast changed abruptly 
in a particular direction, typically by the presence of 
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line segments or edges, but also corners, junctions or 
other singularities (Fig. 1). If the output of the filters 
were statistically independent then we would expect a 
flat correlation as a function of (Ax, Ay, <p, -0). In po- 
lar coordinates (r, 9, <p, tp), the two problems which we 
address are naturally separated: the scaling proper- 
ties result from studying how the histogram depends 
on r (distance) whereas the geometry does it from 
the dependence of the histogram on 9, ip and ip. 

We studied the number of co-occurring pairs of seg- 
ments as a function of their relative distance for dif- 
ferent geometries (Fig 2a,b,c). The different geomet- 
ric configurations correspond to the different orien- 
tations of the segments and their relative position 
within an image. We first studied the number of 
co-occurrences as a function of distance in the line 
spanned by the orientation of the reference segment, 
averaged across all possible orientations of the ref- 
erence line (Fig 2a). When both segments have the 
same orientation we observe a scale invariant behav- 
ior, indicated by a linear relationship in the log-log 
plot (see methods). Also it can be seen from this plot 
that collinear co-occurrences are more frequent than 
any other configuration. Fig. 2b shows the proba- 
bility of co-occurrences is higher for the vertical ori- 
entation, and scale invariance extends over a broader 
range. 

The scaling properties are qualitatively different 
for segments positioned side-by-side, along a line or- 
thogonal to the orientation of the first segment (Fig. 
2c). Iso-oriented pairs were again the most frequent, 
but their co-occurrence in the orthogonal direction 
to the orientation of the first segment (fig 2c, black 
line) does not appear to be scale invariant. This is 
reflected by the presence of a kink as opposed to a 
straight line (power law) in the log- log plot, indicat- 
ing well-defined scales with different behaviour. 

It is worth comparing the scale of interactions one 
observes by using different kinds of filters. Before 
filtering images, the luminance shows correlations, 
which follows a power law behavior J]l], ^o| . After 
applying a Laplacian filter (equivalent to a center- 
surround operator, which measures non-oriented lo- 
cal contrast), the image is mostly decorrelated (Fig. 
2d, red circles) J24|, |2^|. This is seen in the exponen- 
tial decay of the correlations, and the fact that the 



correlations show similar behavior after a pixel-by- 
pixel shuffling of the image (Fig. 2d, cyan circles). 
The strength and scaling of the correlations across 
the collinear line changes radically when one uses an 
oriented filter. In this example, to make a direct com- 
parison between the various filters, we weighted each 
pair of segments by their energy value (linear cross- 
correlation, instead of applying a threshold as was 
done in the earlier calculations). This calculation 
was done for the vertical reference line orientation, 
which showed long-range correlations (Fig. 2b, black 
circles), over much longer distances than observed 
with the Laplacian filter. Moreover, these correla- 
tions were not present when measured in the shuf- 
fled images (Fig. 2, green circles). It is clear from 
the above analysis that when one uses oriented filters 
one reveals strong correlations extending over large 
distances. The next question is how these correla- 
tions depend on the relative orientation of the line 
elements, and whether these dependencies have any 
underlying geometry. We first calculated the total 
number of co-occurrences as a function of the relative 
difference in orientation. Co-occurrences decreased 
as the relative orientation between the pair of seg- 
ments increased, being maximal when they were iso- 
oriented and minimal when they were perpendicular 
(Fig. 3). 

The next observation concerns spatial structure. 
The probability of finding co-occurring pair of seg- 
ments was not uniform but rather displayed a con- 
sistent geometric structure. If the two segments were 
iso-oriented, their most probable spatial arrangement 
was as part of a common line, the collinear configu- 
ration (Fig. 4a). As the relative difference in ori- 
entation between the two segments increased, two 
effects were observed. The main lobe of the his- 
togram (which in the iso-oriented case extends in the 
collinear direction) rotated and shortened, and a sec- 
ond lobe (where co-occurrences were also maximized) 
appeared at 90° from the first (Fig 4a-e). This effect 
progressed smoothly until the relative orientation of 
the two segments was 90°, where the two lobes were 
arranged in a symmetrical configuration, lying at 45° 
relative to the reference orientation. Thus, pairs of 
oriented segments have significant statistical correla- 
tions in natural scenes, and both the average proba- 
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bility and spatial layout depend strongly upon their 
relative orientation. Remarkably, the structure of the 
correlations followed a very simple geometric rule. A 
natural extension of collinearity to the plane is cocir- 
cularity. While two segments of different orientations 
cannot belong to the same straight line, they may 
still be tangent to the same circle if they are tilted 
at identical but opposite angles to the line joining 
them. Given a pair of segments tilted at angles ip 
and (p respectively, they should lie along two possible 
lines, at angles (ip + ip)/2 or (tp + ip + 7r)/2, in order 
to be cocircular (Fig. 4f). This is the arrangement 
we observed in natural scenes. The measured corre- 
lations, given any relative orientation of edges, were 
maximal when arranged along a common circle. To 
quantify this we calculated the orientation of the axis 
where co-occurrences were maximal. We did that for 
different relative orientations and compared it to the 
value predicted by the cocircularity rule (Fig. 5). 
This is particularly remarkable because the compar- 
ison is not a fit, since the cocircularity rule has no 
free parameters. 

Discussion 

We have shown that there are strong, long-range cor- 
relations between local oriented segments in natural 
scenes, that their scaling properties change for dif- 
ferent geometries, and that their arrangement obeys 
the cocircularity rule. The filters we used for edge 
detection in our images were an oriented version of 
Laplacian-like filters, in that they were local but had 
elongated, rather than circularly symmetric, center- 
surround structures. This change is analogous to 
the difference between filters in the LGN and simple 
cells in the primary visual cortex. Thus, given that 
Laplacian filtering decorrelates natural scenes p^| , 
it was surprising to find the long-range correlations 
and scale invariant behavior of the collinear configu- 
ration. It is important to remark that our measure 
of correlation does not differ only on the type of fil- 
ters used (elongated vs. circular symmetric) but also 
on the fact that we measured the correlations along 
a line containing the pair of segments. Long con- 
tours are part of the output of the Laplacian filters 
and thus the image should show correlations which 



might be hidden when integrating them across an 
area, essentially because a curve has zero area and 
thus the correlations along a curve are not signifi- 
cant when integrated over the two-dimensional field 
of view. The findings of long-range correlations of ori- 
ented elements extends the notion that the output of 
linear local oriented filtering of natural scenes cannot 
be statistically independent (2(| and shows that those 
correlations might be very significant through global 
portions of the visual field for particular geometries. 

The cocircular rule has been used heuristically to 
establish a pattern of interactions between filters in 
computer vision Jj], |2?|, |H, pjl, and psychophysical 
studies suggest that the human visual system utilizes 
a local grouping process, "association field", with a 
similar geometric pattern Q . Our finding provides an 
underlying statistical principle for the establishment 
of form and for the Gestalt idea of good continuation, 
which states that there are preferred linkages endow- 
ing some contours with the property of perceptual 
saliency An important portion of the classical 
Euclidean geometry has been constructed using the 
two simplest planar curves, the line and the circle 
]30[ |; here we show that those are, in the same order, 
the most significant structures in natural scenes. 

We have reported the emergence of robust geomet- 
ric and scaling properties of natural scenes. This 
raises the question of the underlying physical pro- 
cesses that generate these regularities. While our 
work was solely based on statistical analysis, we can 
speculate on the possible constraints imposed by the 
physical world. In a simplyfing view, we can think 
of a natural image as composed by object boundaries 
or contours, and textures. Collineal pairs of segments 
are likely to belong to a common contour; thus, our 
finding of scale invariance for collineal correlations is 
in agreement with the idea that scale-invariance in 
natural images is a consequence of the distribution 
of apparent sizes of objects |2l]]. Parallel segments, 
on the contrary, may be part of a common contour 
as well as a common texture, which would explain 
the two scaling regimes we observed. Cocircularity 
in natural scenes probably arises due to the continu- 
ity and smoothness of object boundaries; when av- 
eraged over objects of vastly different sizes present 
in any natural scene, the most probable arrangement 
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for two edge segments is to lie on the smoothest curve 
joining them, a circular arc. These ideas, however, re- 
quiere an investigation which is beyond the scope of 
this paper. 

The geometry of the pattern of interactions in pri- 
mary visual cortex parallels the interactions of ori- 
ented segments in natural scenes. Long-range inter- 
actions tend to connect iso-oriented segments J!?], |l8| 
and interactions between orthogonal segments, which 
span a short range in natural scenes, may be medi- 
ated by short-range connections spanning singulari- 
ties in the orientation and topographic maps in the 
primary visual cortex The finding of a corre- 

spondence between the interaction characteristics of 
neurons in visual cortex and the regularities of natu- 
ral scenes suggests a possible role for cortical plastic- 
ity early in life, in order for the cortex to assimilate 
and represent these regularities. This plasticity might 
be mediated by Hebbian-like processes, reinforcing 
connections on neurons whose activity coincide, i.e., 
their corresponding stimuli are correlated under nat- 
ural visual stimulation. Such plasticity could extend 
to adulthood to accommodate perceptual learning of 
novel and particular forms p^ j. 

While we find coincidences between the pattern 
of interactions in VI and the distribution of seg- 
ments in natural scenes, the sign of the interac- 
tions plays a crucial role. Reinforcing or facilitation 
of co-occurring stimuli (positive interaction) results 
in Hebbian-like coincidence detectors, while inhibit- 
ing the response results in Barlow-like detectors of 
"suspicious coincidences" which ignore frequent co- 
occurrences J3^| . Interestingly the Hebbian idea and 
the decorrelation hypothesis represent two sides of 
the same coin. From our measurements of the reg- 
ularities in natural scenes, and previous studies on 
the higher order receptive field properties in primary 
visual cortex, it appears that both type of operations 
exist. The response of a cell in VI is typically inhib- 
ited when a second flanking segment is placed out- 
side of its receptive field along an axis orthogonal to 
the receptive field orientation. This interaction is re- 
ferred as side-inhibition, which is strongest when the 
flanking segment has the same orientation as the seg- 
ment inside the receptive field [ |l3| , |l5| , Bjj . In the 
present study we found that iso-orientation is the 



most probable arrangement for side-by-side segments 
in natural scenes, which therefore constitutes an ex- 
ample, in the domain of orientation, of decorrelation 
through inhibition. This inhibition may mediate the 
process of texture discrimination 16, [35|. The 
property of end-inhibition has also been interpreted 
as a mechanism to remove redundancies and achieve 
statistical independence |36| . The finding that re- 
sponses of VI neurons are sparse when presented with 
natural stimuli |?7| and models of normalization of 
neuronal responses in VI tuned to the statistics of 
natural scenes J2(| also support the idea that the in- 
teractions in VI play an important role in decorre- 
lating the output from VI. This is consistent with 
the general idea that one of the important functions 
of early visual processing is to remove redundant in- 
39, ^| and suggests that interactions 
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in VI may continue with the process of decorrelation 
which is achieved by Laplacian j24| and local oriented 
filtering 

But the visual cortex also can act in the opposite 
way, reinforcing the response to the most probable 
configurations. This is seen in the collinear configu- 
ration, which is the one that elicits most facilitation, 
and therefore illustrates how VI can enhance the reg- 
ularities in natural scenes. The fact that those corre- 
lations are significant over the entire visual field and 
are highly structured suggests that this is not a resid- 
ual, or second-order process. The opposing processes 
of enhancement of correlations and decorrelation may 
be mediated by different receptive field properties, 
which can exist within the same cell. The same flank 
can inhibit or facilitate depending on the contrast 
|43| , Q suggesting that VI may be solving different 
computational problems at different contrast ranges 
or different noise-to-signal relationship. The dialec- 
tic behavior of visual cortex shows that the interplay 
between decorrelation (extraction of suspicious coin- 
cidences) and enhancement of particular set of regu- 
larities (identification of form) may be mediated by 
the same population of neurons. While the decor- 
relating process may be required to operate in the 
orientation domain to solve the problem of texture 
segmentation, particular sets of coincidences, which 
are repeated in the statistics, such as the conjunction 
of segments that form contours, need to be enhanced 
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in the process of identification of form. 
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Figure Captions 
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Figure 1. An example of the filtering process 
we applied to an image, a) The original image, b) 
The image after processing with local oriented filters 
|p3| . The maximal orientation was calculated at each 
point. The image was converted to binary by consid- 
ering "oriented" only the pixels that after being fil- 
tered at their maximal orientation exceeded a given 
threshold. In the figure, the maximal orientation is 
shown using a color code. 




Distance (Pixels) 



Figure 2. Scaling behaviors for different geomet- 
rical configurations, a) The number of co-occurrences 
between two segments in the relative positions within 
the line that the orientation of the first segment spans 
is shown for different orientations of the second seg- 
ment. This measure was averaged over all possible 
orientations of the first segment. The collinear con- 
figuration is the most typical case and displays a scale 
invariant behavior as indicated by the linear relation- 
ship in the log- log plot, b) The strength of the cor- 
relations and the degree to which it can be approx- 
imated to a power law are more pronounced for the 
particular case in which the reference line segment 
is vertical, c) The same measure when the two seg- 
ments are at a line 90° apart from the orientation 
of the first segment. In all three cases, black corre- 
sponds to iso-orientation, red to 22.5° respect to the 
first segment, green to 45°, blue to 67.5° and yellow to 
90° . d) Full cross-correlation as a function of distance 
for Laplacian filtering (red circles) , oriented filters in 
the collinear vertical direction (black circles) and for 
both cases after shuffling the images. The Laplacian 
filtered image is de-correlated, as can be seen from the 
fact that it shows the same structure as its shuffled 
version (cyan circles). Collinear configuration shows 
long-range correlations, which follow a power law of 
exponent 0.6 (blue line, y=x -0 6 ) and are not present 
when the image is shuffled (green circles). 
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Relative Orientation 

Figure 3. The number of co-occurring pairs of 
segments as a function of their relative difference in 
orientation (£ = ip — (p). These values where obtained 
after integrating the histogram of co-occurrences in 
space for different angular configurations. Each point 
in the graph (£) corresponds to the average and the 
standard deviation of the 16 different configurations 
obtained by choosing one of the 16 possible values 
for the first orientation (ip) and then setting ip = 
(ip + £) (moduli>6) . 



empty figure path 



Figure 4. Plot of the spatial dependence of the 
histogram of co-occurring pairs for different geomet- 
rical configurations, a) The probability of finding a 
pair of iso-oriented segments as a function of their 
relative position (b), a pair of segments at relative 
orientation of 22.5°, (c) 45°, (d) 67.5° or (e) 90°. f) 
Cocircularity solution for a particular example of two 
segments. The solutions to the problem of cocircular- 
ity are two orthogonal lines, whose main have values 
(ip + <ys)/2 or (ip + (p + 7r)/2. For the example given, 



ip(red segment)^ 20°, ip(b\ue segment)^ 40° and the 
two solutions (green lines) are 30° andl20° (all angles 
from the vertical axis). 



Relative Orientation 



■ Natural Scenes 
Cocircularity 



Figure 5. Quantitative analysis of the spatial 
maps. Orientation of the axis where co-occurring 
pairs of oriented elements of relative orientation (£ = 
ip — ip) are maximized. The axis of maximal proba- 
bility was calculated relative to the orientation of the 
segment in the center (tp). This was done for the 16 
possible orientations of ( and the corresponding val- 
ues of (p = (ip + £)(moduld.6), and we computed for 
each ( the mean and standard error. The solid line 
corresponds to the solution predicted by the cocircu- 
lar rule. 



9 



